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Abstract 



Some aspects of the geometry and the dynamics of generaUzed Chaplygin systems are investigated. 
First, two different but complementary approaches to the construction of the reduced dynamics are 
reviewed: a symplectic approach and an approach based on the theory of affine connections. Both 
are mutuahy compared and further completed. Next, a necessary and sufficient condition is derived 
for the existence of an invariant measure for the reduced dynamics of generalized Chaplygin systems 
of mechanical type. A simple example is then constructed of a generalized Chaplygin system which 
does not verify this condition, thereby answering in the negative a question raised by Koiller. 



1. Introduction 



In the past two decades, the rich developments in the field of Geometric Mechanics have led to a 
considerable progress in the study of the geometrical structure and of the dynamics of mechan- 
ical systems with nonholonomic constraints. Several concepts and techniques, familiar from the 
geometric approach to Lagrangian and Hamiltonian mechanics, have been succesfully adapted or 
extended to the framework of nonholonomic mechanics. Examples of this succesful "transfer of 
ideas" can be found, among others, in the symplectic and Poisson descriptions of constrained sys- 
tems |5|, |l5|, ^2], 26, 27, 39 1 and, in particular, in the study of nonholonomic systems with 



symmetry (reduction and reconstruction of the dynamics, stability of relative equilibria ...) where 
elements are being used from the theory of symplectic reduction and from the theory of principal 
connections and Ehresmann connections [^, ^, 2C , |3^ . Whereas most of the foregoing treat- 



ments are concerned with the autonomous case (i.e. time-invariant systems with time-independent 
constraints), the extension to time-dependent nonholonomic systems, using the geometric formal- 



ism of jet bundle theory, has also been treated by several authors: see e.g. |24, 32, 37| and references 
therein. 

An important topic which is receiving growing attention in the literature, concerns the identification 
and characterization of a suitable notion of complete integrability of nonholonomic systems (see 
e.g. j2|, ^, ^, 16, 17, 23, 41]). As is well known, an (unconstrained) Hamiltonian system on a 



2n-dimensional phase space is called completely integrable if it admits n independent integrals 
of motion in involution. It then follows from the Arnol'd-Liouville theorem that, when assuming 
compactness of the common level sets of these first integrals, the motion in phase space is quasi- 
periodic and consists of a winding on n-dimensional invariant tori (see e.g. Q, Chapter 4). For 
the integrability of a nonholonomic system with k constraints one needs, in general, 2n — A; — 1 
independent first integrals. It turns out, however, that for a nonholonomic system which admits an 
invariant measure, "only" 2n — k — 2 first integrals are needed in order to reduce its integration to 
quadratures, and in such a case - again assuming compactness of the common level sets of the first 
integrals - the phase space trajectories of the system live on 2-dimensional invariant tori Q. Several 
authors have studied the problem of the existence of invariant measures for some special classes of 



nonholonomic systems. For instance, Veselov and Veselova [41] have studied nonholonomic geodesic 
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flows on Lie groups with a left-invariant metric and a right-invariant nonholonomic distribution (the 
so-called LR systems). Kozlov [23| has treated the analogous problem for left-invariant constraints. 
Their results have been very useful for finding new examples of completely integrable nonholonomic 
dynamical systems |]l6|, ^ . 

In the present paper we will mainly be concerned with a particular, but important class of nonholo- 
nomic systems: the so-called generalized Chaplygin systems. A system of generalized Chaplygin 
type is described by a principal connection on a principal fiber bundle Q — > Q/G, whose horizontal 
distribution determines the constraint submanifold, and a G- invariant Lagrangian L : TQ — > M. 
Generalized Chaplygin systems exhibit a very nice geometric structure, which has been discussed 
extensively in the recent literature (see e.g. |7, 11, 13, 20,|27, |3^, 43 1, and references therein). One of 
the peculiarities of these systems is that, after reduction, they take on the form of an unconstrained 
system, subject to an "external" force of a special type. In it is shown that in the abelian 
Chaplygin case, the reduced equations can be rewritten in the form of a Hamiltonian system with 
respect to an almost symplectic structure. We will see that this also holds in the nonabelian case. 

In his pioneering paper on the reduction of nonholonomic systems with symmetry, Koiller touches 
the problem of the existence of an invariant measure for the reduced dynamics of generalized 
Chaplygin systems (see ||2^, Section 9). Based on several known examples of such systems which 
do admit an invariant measure, Koiller asks the question whether this property might perhaps 
hold in general. The main result of the present paper is the derivation of a necessary and sufficient 
condition for the existence of an invariant measure for the reduced dynamics of a generalized 
Chaplygin system whose Lagrangian is of pure kinetic energy type. This condition then enables us 
to give a negative answer to Koiller's question by constructing a simple counter example. 

Clearly, there is still much work to be done in order to obtain a deeper and more complete un- 
derstanding of the structure and behaviour of completely integrable nonholonomic systems. Any 
further progress into this matter will also be of importance for the study of perturbations of inte- 
grable nonholonomic systems. Concerning the particular case of generalized Chaplygin systems, it 
would be of interest, for instance, to identify general classes of systems which do admit an invariant 
measure, and to use their characterization to tackle the problem of complete integrability. 

This paper is organized as follows. In Section 2, we first give a very brief review of some aspects 
of the theory of connections on principal bundles before introducing the notion of generalized 
Chaplygin system. Sections 3 and 4 are devoted to two alternative approaches to the description 
of generalized Chaplygin systems, with special emphasis on the reduction and reconstruction of 
their dynamics: a symplectic approach and an approach based on the use of affine connections. In 
Section 5, two examples are treated which exhibit a different behaviour in the reduction process. 
Section 6 then deals with the relation between the symplectic and the affine connection approach. 
The existence of an invariant measure for the reduced equations of generalized Chaplygin systems 
of mechanical type is studied in Section 7. The main result of the paper (Theorem 7.5) presents a 
necessary and sufficient condition for the existence of such an invariant measure, and an example 
is given of a system for which this condition is not satisfied. Finally, some particular results 
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concerning the theory of affine connections are discussed in an Appendix. 

Throughout this paper, we are working in the category of C°°-manifolds (with smooth maps, tensor 
fields, etc.). For convenience, we shah usually not make a notational distinction between a vector 
bundle over a manifold and the module of its smooth sections, i.e. if tt : F ^ N denotes a vector 
bundle over a manifold N (for instance a subbundle of TN), then X £ F simply means that X 
is a smooth section of n. The sole exception will be the sporadic use of the notation X{N) for 
the module of vector fields on N. If V is distribution on a manifold A'^, then its annihilator is a 
codistribution on denoted by T>". Both T> and will also be identified with the corresponding 
vector subbundles of TN, resp. T*N. Finally, the tangent map of a mapping / between manifolds 
will be denoted interchangeably by Tf and f^. 

2. Generalized Chaplygin systems 

In this section, we first briefly review the notion of a principal connection on a principal fibre bundle. 
For details we refer to ||l^ . Secondly, we consider a geometric framework for general nonholonomic 
mechanical systems and then introduce the concept of generalized Chaplygin systems. 

2.1. Principal connections 

Let Q be the configuration manifold of a physical system and assume that there is a left action of 
a Lie group G on Q: 

^ : GxQ — > Q 

{g,q) ^ q,{g,q) = q,g{q) = gq. 

Note that we consider here left actions, which is the usual convention in mechanics, instead of the 
right ones, as considered for instance in [^]. The orbit through a point q is OrbG('?) = {gq \ g S G}. 
We denote by Q the Lie algebra of G. For any element ^ € 0, will denote the corresponding 
infinitesimal generator of the group action on Q. Then, 

T,iOThGiq)) = {CQiq)\CeQ}. 

Assuming that the action ^' is free and proper, we can endow the quotient space Q/G = M with 
a manifold structure such that the canonical projection vr : Q — > M is a surjective submersion. 
In the framework of the mechanics of (coupled) rigid bodies, for instance, the quotient manifold 
M is commonly called the shape space of the system under consideration. We then have that 
Q{M,G,tt) is a principal bundle with bundle space Q, base space M, structure group G and 
projection vr. Note that the kernel of •7r*(= Tit) consists of the vertical tangent vectors, i.e. the 
vectors tangent to the orbits of G in Q. We shall denote the bundle of vertical vectors by Vn, with 
(v.), = r,(OrbG(g)), q G Q. 

A principal connection on Q{M, G, vr) can be defined as a distribution TC on Q satisfying: 
(i) TgQ = e (V.)„ yq G Q; 
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(ii) Hgq = Tq"^g{TCq), i.e. the distribution H is G-invariant; 

(iii) Hq depends smoothly on q. 

The subspace TCq of TqQ is called the horizontal subspace at q determined by the connection. 
Alternatively, a principal connection can be characterized by a g-valued 1-form 7 on Q satisfying 
the following conditions 

(i) li^Qiq)) = ^ for ah e G 0, 

(ii) j{T^gX) = Adg{j{X)) for ah X G TQ. 

The horizontal subspace at q is then given by Hq = {vq £ TqQ \ 7(^5) = 0}. A vector field X on 
Q is called horizontal if X{q) £ Jiq at each point q. 

Given a principal connection, we have that every vector v G TqQ can be uniquely written as 

V = Vl + V2 , 

with vi G TLq and V2 G (Vjrjg. We denote by h : TQ — > Ti and v : TQ — > Vtt the corresponding 
horizontal and vertical projector, respectively. The horizontal lift of a vector field y on M is the 
unique vector field on Q which is horizontal and projects onto Y . 

The curvature 17 of the principal connection is the 0- valued 2-form on Q defined as follows: for 
each q £ Q and u,v £ TqQ 

n{u, v) = d7(hn, hv) = -7([C/^ V\) , 

where and are the horizontal lifts of any two (local) vector fields U and V on M for which 
U^{q) = h.u and V^{q) = hf, respectively. The curvature measures the lack of integrability of 
the horizontal distribution and plays a fundamental role in the theory of holonomy (see ||l^ for a 
comprehensive treatment). 

2.2. Nonholonomic mechanics and generalized Chaplygin systems 

A nonholonomic Lagrangian system consists of a Lagrangian L defined on TQ, the tangent bundle 
of an n-dimensional configuration manifold Q, and constraints which determine a submanifold T> 
of TQ. This means that the only allowable velocities are those belonging to T>. In case P is a 
vector subbundle of TQ, we are dealing with the case of linear constraints and if, in addition, 
this subbundle corresponds to an integrable distribution, we are reduced to the case of holonomic 
constraints. In the sequel, we will always assume that tq{V) = Q, where tq : TQ — > Q is the 
tangent bundle projection. 

Let (q^), A = 1, n, be local coordinates on Q and denote the induced bundle coordinates on TQ 
by {q^, q^). In a local description, a constraint submanifold T> of codimension k can be defined by 
the vanishing of k independent functions (pi (the constraint functions). In the following, we shall 
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only consider the case of linear constraints, such that the functions 0j can be taken to be of the form 
4>i{(ljQ) = MiAC?)^"^, ^ = 1, ■■■,k. Application of d'Alembert's principle, using the common notion 
of virtual displacement for the case of linear constraints, then leads to the constrained equations 
of motion 

which, together with the constraint equations fJ,iAQ^ = {i = 1, determine the dynamics of 

the nonholonomic system. Here, the A' are Lagrange multipliers to be determined. The right-hand 
side of equation (||) precisely represents the "reaction force" induced by the constraints. 

We now describe the structure of a so-called generalized Chaplygin system [0, |2^, The 
configuration manifold Q of a generalized Chaplygin system is a principal G-bundle vr : Q — > Q/G, 
and the constraint submanifold V is given by the horizontal distribution Ti. of a principal connection 
7 on vr. Furthermore, the system is described by a regular Lagrangian L : TQ — M, which is 
G-invariant for the lifted action of G on TQ. In this paper, we shall mainly restrict our attention 
to systems of mechanical type for which L = T — V , where T : TQ — > M is the kinetic energy, 
corresponding to a Riemannian metric g on Q, and V : Q — > R is the potential energy. Whenever 
we add the word "mechanical" to the description of a system, we shall always be referring to this 
situation. For the case of a Chaplygin system we then suppose, in addition, that both the potential 
energy and the metric g are G-invariant so that 

for all ^ G 0. In particular, it follows that all fundamental vector fields are Killing vector fields. 

Typical problems in mechanics, such as the vertical and the inclined rolling disk, the nonholonomic 
free particle and the two wheeled carriage, can be interpreted as generalized Chaplygin systems in 
the above sense. Systems of that type also occur in many problems of robotic locomotion and 
motions of microorganisms at low Reynolds number |3^] . The dynamics of generalized Chaplygin 
systems can be described from a symplectic point of view, as we will show in the next section. 
However, for generalized Chaplygin systems of mechanical type there also exists a nice geometric 
description in terms of affine connections. This will be outlined in Section 4. 

We conclude this section with some comments concerning terminology. Classically, a mechanical 
system with Lagrangian L{q^,q^), A = l,...,n, subject to k linear nonholonomic constraints, 
is said to be of Chaplygin type if coordinates {q'^,q") can be found, with a = 1,. . . ,k and a = 
k + 1, . . . ,n, such that the constraints can be written in the form q"' = B^{q^~^^, . . . , q^)q°^ and 



such that L does not depend on the coordinates q°' (see e.g. |33]). Such a system can be (locally) 



interpreted as a special case of the generalized Chaplygin systems introduced above, with Q = M"' 
and with an action defined by the abelian group G = W (cf. @). Koiller refers to the more general 
case, considered in the present paper, as "non-abelian Chaplygin systems". In the literature on 
nonholonomic systems with symmetry, these systems are also said to be of "principal" or "purely 
kinematical" type |Q, 10, 20 1. Finally, it should be emphasized that the generalized Chaplygin 
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systems studied in [ 27 [ are still of a more general type than the ones we consider here in that they 
are defined on fibre bundles which need not be principal bundles. 



3. Symplectic approach 



As mentioned in the Introduction, considerable efforts have been made to adapt and extend several 
ideas and techniques from the geometric treatment of unconstrained problems to the study of 
systems with nonholonomic constraints. The subject has been approached from different points 
of view: a Lagrangian approach a Hamiltonian approach |l^] and a formulation in 



terms of (almost-)Poisson structures |12, 22, 30, 39|. In this section, we start with a brief review 



of an interesting symplectic approach to nonholonomic dynamics, developed in ^] which, in 
particular, is well suited for the treatment of nonholonomic systems with symmetry. We first 
outline this approach for general nonholonomic systems, and then turn to the case of generalized 
Chaplygin systems for the discussion of reduction and reconstruction of the constrained dynamics. 

We start by fixing some notations. In terms of the tangent bundle coordinates {q^,q'^), let us 

d d 
denote hy A = the dilation vector field on TQ and by S" = dq^ (8) 7—^ the canonical vertical 

oq^ oq^ 

endomorphism (see |2^). The action of S on a 1-form will be denoted by S*. Then we can define 
the Poincare-Cartan 1-form and 2-form, corresponding to a given Lagrangian L, by 9^ = S*dL 
and lvl = —dOi, respectively. We further have that Ej^ = AL — L represents the energy function 
of the system. If the Lagrangian L is regular, which will always be tacitly assumed in the sequel, 
lol is symplectic and induces two isomorphisms of C°°(TQ)-modules: 

where bL(X) = ix^^L and (II = ^l^- In the absence of constraints, the dynamics of the Lagrangian 
system, with Lagrangian L, is given by the (unique) solution Tl of the equation ir^^L = dE^, i.e. 
Tl = thidEi). Indeed, Tl is a second-order differential equation field (SODE for short) whose 
integral curves {q^{t).,q^{t) = ^^%-{t)) are determined by the solutions q^{t) of the Euler-Lagrange 
equations for L. 

In the presence of nonholonomic constraints, the equations of motion must be modified in order 
to incorporate the constraints into the picture. Since we confine ourselves here to the case of 
linear constraints, the constraint submanifold P is a vector subbundle of TQ, determined by a 
(nonintegrable) distribution on Q which we also denote by T>. In addition, we shall always assume 



that the constraints verify the so-called "admissibility condition" (see e.g. [27|), i.e. for all x € 

dim(r,P)° = dim5*((r,P)"), 
where the annihilator of TxT) is taken in T*TQ. Locally, D is described by equations of the form 

fJ'iAiq)q^ = 0, i = l,...k, 
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with k = codim (2?). Next, we define a distribution F on TQ along V, by prescribing its annihilator 
to be a subbundle of T*TQ\p which represents the bundle of reaction forces. More precisely, we 
set F° = S*{{TT>)°). The equations of motion for the nonholonomic system are then given by 

f {ixOJL - dEL)\v G F° , 
\ G TP . 

This system will have a unique solution X provided the "compatibility condition" holds, i.e. fl 
TT> = 0, where F-^ = ^l{F°). For systems of mechanical type, with linear nonholonomic constraints 
and positive definite kinetic energy, this condition is always fulfilled. Indeed, putting 

dj = -jJ-iAW^^fijB, i,j = l,...k, 

where (W^^) is the inverse of the Hessian matrix ( g -aq -b j ' compatibility locally translates 
into regularity of the matrix {Cij) (cf. |2^, where the compatibility condition for a nonholonomic 
system was called the regularity condition). Under this condition, one can show that we have a 
direct sum decomposition of Tx>TQ = TV © F-^ and that the constrained dynamics X is obtained 
by projecting the unconstrained Euler-Lagrange vector field Tl (restricted to V) onto TV with 
respect to this decomposition. 

It should be pointed out that the solution X of (P) satisfies automatically the SODE condition 
along V, i.e. S{X)^d = Ap. This implies that, in local coordinates, the integral curves of X 
on V are of the form {q"^ (t) , q"^ (t) = ^^(t)), whereby the q^{t) are solutions of the system of 
differential equations (|^), together with the constraint equations fj.iA{q)q'^ = 0, i = The 
local coordinate expression for X reads 

X = q^ ^ 



dq^ 

where (C*-^) is the inverse of the matrix {Cij) introduced above and where, for ease of writing, we 

have put pA = Trrj^ ^ = 1) •••>'^- 

oq^ 

Before proceeding, we recall that for nonholonomic Lagrangian systems with constraints which are 
linear (or, more general, homogeneous) in the velocities, the energy El is a conserved quantity (see 
e.g. [|lO|). This, therefore, in particular applies to the generalized Chaplygin systems considered in 
this paper. 

3.1. Reduction 

The reduction theory of nonholonomic systems with symmetry, and related aspects, has become 
an intensive field of research (5|, 0, |ll|, 13, 22, 30, 34]. Here we consider the special case of 
nonholonomic systems of generalized Chaplygin type. 

The given data are (cf. Section 2): a principal G-bundle vr : Q — > M = Q/G, associated to a free 
and proper action ^ of G on Q, a Lagrangian L : TQ — > M which is G-invariant with respect to the 
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lifted action on TQ, and linear nonholonomic constraints determined by the horizontal distribution 
(here denoted as 2?) of a principal connection 7 on tt. Taking into account the available symmetries, 
we can reduce the number of degrees of freedom of the problem. In the following, we review the 
various geometric concepts involved in the symplectic approach to this reduction process |ll|, |l^ . 

The symplectic action. Consider the lifted action of G on TQ, i.e. : G x TQ — > TQ with 
^{g,Vq) = T'$g{vq){= "^g^i^Vq)) for any g £ G and Vq G TqQ. This action is also free and proper 
and, moreover, it is symplectic with respect to ujl. For any ^ G 0, the infinitesimal generators ^tq 
and of ^ and respectively, are rg-related, i.e. 

TQ* ° Ctq = ^q°tq . (3) 

Let us denote by p : TQ — > TQ = TQ/G the natural projection. From the given assumptions 
it follows that the constraint submanifold T>, the energy and the vector subbundle F are G- 
invariant. The induced action of G on V, i.e. the restriction oi ^ to G x T>, is still free and proper 
and we can regard the orbit space V = T>/G as a submanifold of TQ/G. Note that there exists a 
natural identification V = Q Xq/qT{Q/G) as principal G-bundles over T{Q/G). The isomorphism 
is obtained by mapping Vq (z T> onto (g, vr^, (vq)). It then follows that V/G can be naturally identified 
with T{Q/G) and we have, in particular, 

P\v = '^*\v- (4) 

Henceforth, the restriction of /? to P will also be simply denoted by p. 

The connection. A direct computation shows that along D we have that VpCi F 
denotes the subbundle of TTQ which is vertical with respect to the projection p, i.e 
Observing that Vp|p C TV, it then easily follows that 

TV={Fn TV) e Vpp . 

Since U = F f] TV is G-invariant, the above decomposition of TV defines a principal connection 
T on the principal G-bundle p : V — > V = V/G, with horizontal subspace Ux at each x (z V. 
Moreover, we have that the dynamics X of the generalized Chaplygin system, determined by @, 
belongs to U. We shall denote the projection of the bundle U under Tp{= p*) by U. 

We can summarize the situation in the following diagram: 

TV = U®Vp ^ TV^U 

i i 

V T{Q/G) ^ V 

The principal connection T is obviously related to the original connection 7 of the Chaplygin 
system. Indeed, take w G T^^V and consider tq^w G TqQ. Then, we can write 

TQ*w = {tq^w - i7iTQ^w))Q{q)) + i7iTQ^w))Q{q) , 



= 0, where Vp 
. Vp = KerTp. 
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where tq^w - {-f{TQ^w))Q{q) G Vg and {7{TQ^w))Q{q) G (V^)g. Putting -firQ^w) = ^ G g, a direct 
computation shows that w — CxQivq) G ^ and, consequently, w = {w — S,TQivq)) + CTQ{vq) is the 
{U,Vp) decomposition of w. Herewith we have proved the fohowing property. 

Proposition 3.1. The connection 1-forms T and 7 are related by T = Tq'j, i.e. T^^{w) = 
lq{TQ*w) for any Vq e V and w G T^^V. 

Let us denote the horizontal projectors, associated to 7, resp. T, by : TQ — > V{c TQ), resp. 
hY : TV — > U{c TV). Likewise, the vertical projectors onto Vtt, resp. Vp, will be denoted by v^, 
resp. vx. In order not to further overload the notations, we will use the same superscript ^ for the 
horizontal lifts of vectors (vector fields) with respect to either 7 or T; in principle it should always 
be clear from the context which horizontal lift operation is being used. 

We now have that 

i.e. the following diagram is commutative: 

{TTQ d) TV ^ TV (C TTQ) 

TQ ^ TQ 

Indeed, taking into account Proposition 3.1 we see that for any w G T^^V, 7g(rQ^hx(w)) = 
'^vqiMxiw)) = and, hence, rg^h-f (if) is horizontal with respect to 7, i.e. TQ^h-f{w) G Vg. By 
definition we also have h^{TQ^w) G Vq. Using the fact that vr o tq = tqjq o tt^ we obtain: 

7r*(h^(TQ^u;)) = Ti^TQ^w = Tq/g^{'k*)*w = Tq/g^p^w, 

where the last equation follows from Similarly, we have: 

We thus see that the 7-horizontal tangent vectors at q, tq^ (h-fw) and h^(rQ^tt;), have the same 
projection under vr* and, therefore, they are equal. This completes the proof of (|5|). Denoting the 
curvature tensors of the principal connections 7 and T by W and Vt^ , respectively, one can easily 
deduce from Proposition 3.1 and (^) the relation 

= t*qW. (6) 

The 1-form. Denote by 6' the pullback to V of the Poincare-Cartan 1-form 9^, i.e. 9' = j^9L, 
where j-p :V ^ TQ is the canonical inclusion. By means of the solution X of (||) we can construct 
a 1-form ax on V as follows: 

ax = ix {hr*d9' - dhr*9') , (7) 

with the usual convention that, for an arbitrary p-form /? on V, h-f*/? is the p-form defined by 
hr*/?(Xi, ...,Xp)= /3(h-f (Xi), . . .,hr{Xp)). 
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The Lagrangian. The Lagrangian L of the given mechanical system induces a Lagrangian L* : 
T{Q/G) — > R on the quotient space V ^ T{Q/G), given by L*{q,Vq) = L{q,v^) for any q € 
'K^^{q), and where denotes the 7-horizontal hft of Vq at q. This is well-defined because of the 



G-invariance of L. Moreover, one can show that L* is a regular Lagrangian on T{Q/G) (cf. |27|). 



Now, we are in a position to state the following reduction result (see |11, pO[) 



Proposition 3.2. The dynamics X of the generalized Chaplygin system projects onto T>, and its 
projection X is determined by the equation 

ix^L* = dEi* + ax , (8) 
where 'ax is the projection of the 1-form ax, defined by (j^. Moreover, we have that ixOx = 0. 

Remark 3.3. The sign of ax in equation (^) differs from the one in the corresponding expression 
derived in [11|, where the discussion took place in a more general symplectic framework with an 
exact symplectic structure uj = dO (whereas here we have ujl = —dOi). The signs would agree if 
we would have defined 9' as —j^Oi. 

It can be easily verified that the form ox is a semi-basic 1-form on P (see also Section 6), from 
which it then follows that the vector field X, defined by (^), is a SODE. Moreover, one can show 
that not only the contraction of ox with X vanishes, but that, more generally, iyox = for any 
SODE Y on "D. We thus see that a generalized Chaplygin system reduces to an unconstrained 
mechanical system, with an external nonconservative force of "gyroscopic" type, which is geomet- 
rically represented by the 1-form ZSx (see also pO| , [3^ ) . The "gyroscopic" character of this force is 
also in agreement with the fact that the projected energy function E^* is a conserved quantity of 
the reduced dynamics. But there is more to be said about it. 

Proposition 3.4. The 2-form S = \v^*dO' - dh.-r*6' on V projects onto a 2-form T. onV and the 
1-form Hx satisfies 

OLx = ■ixS. 

Proof. Let be the fundamental vector field of the G- action on D, induced by an arbitrary 
element G 0. We must prove that ^x> belongs to the characteristic distribution of the 2-form S. 
First, we have that 

= —i^^dYi-f* 0' . 

For any vector field Y on V, i^-^dhr* 9' {Y) = iv{9' {"tirY)) - 9' {\ir[iv ,Y]) . Now, if Y is vertical, 
we readily see that i^^dli^* 9' iX) = ^- If ^ is horizontal, we have 

i^^dhr*9'{Y) = C^J'{Y) = Q, 

because of the G-invariance of 9' . It therefore remains to prove that i^^dT, = 0. For any two vector 
fields Y, Z we have that 

i^^d^{Y,Z) = {i^^dhr*d9'){Y,Z) 

= {hr*d9') {Y, Z) - \ir*d9'{[iv, YlZ) + hr*d9'(\iv, Z],Y) . 
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If at least one of the vector fields Y and Z is vertical, then i^^dJ^{Y, Z) = 0. Taking Y and Z both 
horizontal, we find, taking into account the G-invariance of dO' and (,x> £ ^(^)) 

Z) = iv {de'{Y, Z)) - de'{[Cv, Y],Z)- de'{Y, [Cv, Z]) = {C^^de'){Y, Z) = . 

The last part of the proposition now immediately follows from (^) and the projectability of X. 

\qed\ 

Consider a local trivialization ?7 x G of vr, with coordinates where i = I, . . .k = dimG 

and a = 1, . . . ,n — k. Choosing a basis e, (i = 1, . . . , k) of the Lie algebra Q, and using the left 
trivialization TG = G x g, a tangent vector v € T(^ g)(C/ x G) = T^U x g can be represented by 
a pair {w,S,), whereby w £ T^U and ^ = ^*ej G g. In terms of the coordinates {q^ , , q"" , C) on 
T{U X G) the G-invariant Lagrangian can then be written as 

Strictly speaking, I represents the reduction of L to TQ/G. With respect to the given local 
trivialization, we further denote the connection coefficients of the given principal connection 7 by 
= rjj(g^, . . . , q"'~^), and the constraints then take the form ^* = —T\q°'. In particular, it follows 
that the reduced Lagrangian L* is given by L*{q°-,q"-) = £{q"',q°', —Tlq^). 

With all the above one can now derive the following coordinate expression for the reduced dynamics 
(see also 0, |o|): 



(d^L* \ dL* 
— TT-r 1 1 Pa = -T^ — 1 and audq^ is the local 
dq^dq^ J dq" 

expression for the gyroscopic 1-form ax . The ah are explicitly given by 



(^eO {dq^ d<f 

where the * on the right-hand side indicates that, after computing the derivative of H. with respect 
to one replaces the everywhere by —Taq"'- The constants c*-^ appearing in the last term on the 
right-hand side are the structure constants of g with respect to the chosen basis, i.e. [cj, e^] = c*^ej. 

Note in passing that the expressions -^-^ — -^-^ — c*,^.r;^r^ are the coefficients of the curvature of 
7 in local form. 

3.2. Reconstruction 

A natural problem related to the reduction of mechanical systems with symmetry concerns the 
reverse procedure: once the solutions of the reduced dynamics have been obtained, how can one 
recover from it the solutions of the original system. This is called the "reconstruction problem" of 
the dynamics. This problem is intimately related to the concepts of geometric and dynamic phase, 
which play an important role in various aspects of mechanics pl| and in the study of locomotion 



systems (for example, in the generation of net motion by cyclic changes in shape space |18, 35 1) 
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In the case of a generalized Chaplygin system, the reconstruction of the dynamics on V from the 



reduced dynamics on V proceeds as fohows |13]. Suppose that the flow of the reduced system X 
is known. Let c{t) be the integral curve of X starting at a given point x € P, and fix j; S p~^{x). 
We want to find the integral curve c(t) of X with initial point x and projecting onto c{t), i.e. 
p{c{t)) = c{t). But this is precisely the horizontal lift through x of c(t), with respect to the 
principal connection T. We recall here the proof of this simple fact (see also jl3[|). 



Proposition 3.5. The integral curve c{t) of the generalized Chaplygin system, starting at x V, 
is the horizontal lift with respect to the principal connection T of the integral curve c{t) of the 
reduced system starting at x = p{x) . 

Proof. Let d{t) denote the horizontal lift of c(t), starting at x (say, at t = 0). Therefore, d{t) G 
T) and p{d{t)) = c{t) for each t, and (i(0) = x. Since X and X are p-related, we have that 
Tp{X{d{t))) = X{p{d{t))) = X{c{t)) = Tp{d{t)). Hence, d{t) - X{d{t)) is vertical for each t. On 
the other hand, it is also horizontal, since X gU. Therefore, we deduce that d{t) = X{d{t)). 

\qed\ 

Summarizing: in the symplectic formalism, the reconstruction simply consists of a horizontal lift 
operation with respect to the induced connection T living on V. 



4. Affine connection approach 



In this section we will describe the reduction and reconstruction problem for generalized Chaplygin 
systems of mechanical type, from a different point of view, namely in terms of the theory of affine 
connections. First, we recall some general results which apply to any kind of nonholonomic system 
of mechanical type. 

Let Q be the configuration space of a mechanical system with Lagrangian of the form 

Hv) = ^g{v, v) - V o tq{v), V G T^Q , 

where g' is a Riemannian metric on Q and V : Q — > M is the potential energy function. We denote 
by the (covariant derivative operator of the) Levi-Civita connection associated to the metric 
g. It is well-known that a curve c : I — > Q is a solution of the Euler-Lagrange equations for the 
(unconstrained) Lagrangian if and only if 

V^(^)C(t) = -grady(c(t)), 

where the gradient is also considered with respect to the metric g. 

For nonholonomic systems there is a similar description. Let V again denote a (nonintegrable) 
distribution on Q, describing some linear nonholonomic constraints. The second-order differential 
equations (jl]) for the mechanical system with Lagrangian L and constraints V, can be written 
intrinsically as 

V^(t)C(i) + grady(c(t)) G P,^,) , c{t) G , 
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where here denotes the ^-orthogonal complement to V (see e.g. |2^, 38, 4C]). 
Alternatively, if we denote by 

V:TQ^V, Q :TQ — > , 
the complementary g'-orthogonal projectors, we can define an affine connection 

= v^y + (v^Q)(y), 

such that the nonholonomic equations of motion can be rewritten as 

Vc(t)c(t) = -V{giadV{c{t))), 
and where we select the initial velocity in V (cf. [p^]). 

In what follows, we shall restrict our attention to Lagrangians of "pure kinetic energy type", i.e. 
we assume V = 0. The reason for doing this is twofold. First, it makes the geometric picture more 
clear and tractable, in that the equations of motion for the nonholonomic mechanical system can 
then be seen as the geodesic equations of an affine connection. Secondly, the extention to systems 
with a nontrivial potential energy function is rather straightforward but, at least for those aspects 
of nonholonomic dynamics that are of interest to us here, it does not really tell us anything new. 

It can be easily deduced from its definition that the connection V restricts to P, that is, 

for all y € 15 and X £ X{Q). This kind of affine connections, which restrict to a given distribution, 
have been studied in In particular, such a behaviour implies that the distribution V is 

geodesically invariant, that is, for every geodesic c{t) of V starting from a point in D, c(0) € 
^c(o)) have that c{t) € I5c(t)- [H]) & iiice property is derived which characterizes geodesically 
invariant distributions in terms of the so-called symmetric product of vector fields, which is defined 
as {X : Y) = + VyX. This property asserts that P is geodesically invariant if and only if we 

have that {X -.Y) £ D, VX, Y £ V. Note in passing that the symmetric product of vector fields is 
a differential geometric concept with important applications to control theory, first introduced in 



|29|. For instance, it plays a fundamental role in the controllability analysis of mechanical control 



systems ||9|, |T^, 29], in the description of the evolution of these systems when starting from rest 
and in the design of motion planning control algorithms 

Proposition 4.1. For all Z € X{Q) and X , Y e V we have that 

Z {g{X, Y)) = g{VzX, Y) + g{X, V zY) . 

Proof. In view of the definition of V, we have that 

g{VzX, Y) + g{X, V zY) = g{V%X, Y) + g{X, V|y) + <7((V|Q)(X), Y) + g{X, (V|Q)(y)) 

= Z{g{X,Y)) , 
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since (V|Q)(X), (V|Q)(y) G (see Proposition 6.1 in fH). \qed\ 

We derive from this proposition that the connection V also has the following property: parallel 
transport is an isometry along the distribution T). 

A direct computation shows that the torsion of V is the skew-symmetric (l,2)-tensor field 

f{X,Y) = {V%Q){Y)-{V^yQ){X). 

Observe that if X, y G D then f{X, Y) ^V^. 
4.1. Reduction 

We now return to the case of generalized Chaplygin systems of mechanical type, but again assuming, 
for simplicity, that L only consists of a kinetic energy part (i.e. V = Q). From the above we then 
know that the equations of motion of the system can be written as 

V,(i)c(t) = 0, c(0) G V. 

The following is mainly inspired on the reduction process of generalized (or non-abelian) Chaplygin 
systems as described in [pO| . 

Let us define a metric g on the base manifold M(= Q/G) as follows 

gx{Ux,Vx) = gq{Uq,Vq), X G M, Ux,Vx GT^M , 

where q G 7r~^{x) and Ug, Vq are horizontal vectors which project under vr onto Ux and Vx, respec- 
tively. From the G-invariance of g we deduce that the right-hand side is independent of the chosen 
point q in the fibre 7r~^(x) and, hence, g is well defined. 

Proposition 4.2. We have that for all X,Y £ X(M) and ^ G 

'Csg(V^.y")=0. 



Proof. Since .fg is a Killing vector field, it follows from Proposition VI. 2.2 in ||19[ that 

v^,y^ = ,,,jy^ + v^, [^q, y'^] = o , (lo) 

because Y'' and X'' are projectable. Therefore, we only need to prove that 



This condition is equivalent to 



0. 



Q(v^,y'^) =0. (11) 



But, as "D and P are G-invariant, we have that C-^^Q = 0. This, together with (|10|), imply (pTf). 

\QEd\ 

Now, we define an affine connection on M as follows: for X, y G X(M), put 

v^y = vr*(Vx'.y'^). 
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This is well-defined since, by Proposition 4.2, the vector field Vxh^'* is projectable, and one easily 
verifies that V satisfies the properties of an affine connection. Then, we obtain the following 
important result. 



Proposition 4.3.(|^0[) The geodesies ofV, with initial condition in V, project onto the geodesies 
ofV. 

Proof. Key fact for the proof is that 2? is geodesically invariant with respect to V. \QED\ 

Consequently, we have found that the equations of motion of the given generalized Chaplygin 
system reduce to the geodesic equations of the induced affine connection V on M = Q/G. 

Consider the following (0,3)-tensor field on Q: 

where is the horizontal projector and ^l"^ is the curvature of the connection 7. Observe that K 
is horizontal, i.e. it vanishes if one of its arguments is a vertical vector, and it is skew-symmetric 
in its last two arguments. Moreover, one can see that 

RggiT-i^giUg), T^g{Vg), T^g{Wg)) = Kg{Ug, Vg , Wg) , 

for all g (z G and q (z Q. Consequently, K induces a (0,3)-tensor on the base manifold M 

Kx{Ux,Vx,Wx) = Kg{Uq, Vq, Wq) , 

where 7r(g) = x and Uq, Vq, Wq are tangent vectors in q projecting onto Ux,Vx, Wx, respectively. K 
(resp. K) is called the metric connection tensor on Q (resp. M). 



In |2[1[ |, it was shown that application of the so-called Hamel's approach to mechanics, leads to two 
additional affine connections on M, whose geodesies are also solutions of the reduced nonholonomic 
problem. These connections are given by 

{V^)xY = V^^Y + B{X,Y) , {V^)xY = V^xY + B{Y,X), (12) 

where B is the (l,2)-tensor field defined by p{B{X,Y)) = K{X,Y,^g/3), for any /3 G A\Q/G), X, 

Y € X{Q/G). In general, the tensor which encodes the difference between an affine connection 

V on a Riemannian manifold and the Levi-Civita connection corresponding to the Riemannian 
metric, is called the contorsion of V. So, B here represents the contorsion of V{^. 

The following explicit formula for the connection V was then derived in ||2^ (up to a minor 
misprint): 

VxY = V^y + ^ {B{X, Y) + B{Y, X) - C{X, Y)) , 

where G is the (l,2)-tensor field implicitly defined by j3{G{X,Y)) = iir(tJg/?, X, y), for arbitrary 
(5 G A^{Q/G), X,Y e X{Q/G). As noted in [|0|, the average of Hamel's connections, i.e. V^/^ = 
^(V(^ + ^2)1 general differs from V, because of the skew-symmetric term C(X, Y). 
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It is interesting to observe that from Proposition 4.1 one can easily deduce 



Vg = 0, 

that is, V is a metric connection. From Proposition A. 3 (see Appendix) and the definition of 
-B, it is readily seen that is also a metric connection. In general, however, the connections S/2 
and V^/^ will not be metric. In fact, it is not hard to prove the following result. 

Proposition 4.4. The following properties are equivalent: 

(i) V2 metric; 
(11) is metric; 

(Hi) the tensor field B is skew- symmetric; 
(iv) V^/^ is the Levi-Civita connection of g. 

Later we will see that these properties are also equivalent to the vanishing of the 1-form Wx and, 
hence, to the Hamiltonian nature of the reduced system (cf. Corollary 6.2). 

The torsion of V is given by T{X, Y) = vr* T{X^ , Y^), with T the torsion of V. Then, we see from 
the above that the metric connection V is the Levi-Civita connection associated to g iff the torsion 
of V takes values in the vertical tangent bundle to vr for each pair of vectors in V. 

Finally, the following result shows that equality of V and V^^^ is a rather strong condition. 

Proposition 4.5. 

Proof. If V = V''^/^, then V^/^ is metric. By Proposition 4.4, this implies that V^/^ coincides 
with V^. The reverse implication is trivial. \QED\ 

4.2. Reconstruction 

In view of Proposition 4.3 above, we see that, in the present setting, the reconstruction of the 
solution curves in Q of the given constrained system, from those of the reduced system on M, 
consists of a horizontal lift operation with respect to the connection 7. That is, let c{t) be a 
geodesic of V and choose c(0) G Q such that 7r(c(0)) = c(0). Then, the geodesic starting at c(0), 
with initial velocity c(0) S ^^0(0)1 is precisely the horizontal lift of c{t) with respect to the principal 
connection 7. 



5. Examples 



5.1. Mobile robot with fixed orientation 

Consider the motion of a mobile robot whose body maintains a fixed orientation with respect to 
its environment (see |18] for more details). The robot has three wheels, with radius i?, which turn 
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simultaneously about independent axes, and perform a rolling without sliding over a horizontal 
floor. Let (x, y) £ denote the position of the center of mass of the robot (in a cartesian 
reference frame, with horizontal x- and y-axis), 6 £E>^ the steering angle of the wheels and ip €E>^ 
the rotation angle of the wheels in their rolling motion over the floor. The conflguration space can 
then be modelled by Q = §^ x §^ x M^. 

The Lagrangian function L is the kinetic energy function corresponding to the metric 

g = mdx ® dx + mdy ® dy + JdO iS> d9 + 3Jwdip (8) dip , 

where m is the mass of the robot, J its moment of inertia and the axial moment of inertia of 
each wheel. The constraints are induced by the conditions that the wheels roll without sliding, in 
the direction in which they "point" , and that the instantaneous contact points of the wheels with 
the floor have no velocity component orthogonal to that direction (cf. p^]): 

X sin 9 — y cos 9 = , 
X cos 9 + y sin 9 — Rip = . 

The constraint distribution T> is then spanned by 

/A A+i? fcos0— +sin0— H 
\d9 ' dtp \ dx dy J ) 

If we consider the abelian action of G = on Q by translations 

^ : GxQ — > Q 

{{a,b),{9,ip,x,y)) i — > {9,'tp, a + x,b + y) , 

we see that the constraint distribution V can be interpreted as the horizontal subspace of the 
principal connection 7 = (dx — Rcos6d'ip)ei + (dy — Rsm9d'ip)e2, where {61,62} is the canonical 
basis of (identified with the Lie algebra of G). 

The metric induced on M = Q/G here becomes 

g = Jd9 <^d9 + (SJ^; + mR^)dip (g) dip . 

The reduced Lagrangian L* is the kinetic energy function corresponding to the metric g. Moreover, 
one easily verifies that the gyroscopic 1-form ZSY identically vanishes and, hence, the symplectic 
reduction (^) yields 

i.e. the reduced system is an unconstrained, purely Lagrangian system. 

On the other hand, one can easily check that, in this example, the metric connection tensor K also 
vanishes. Consequently, with the notations of the previous section, V = \7^^^ and by Proposition 
4.5 we then have that V = V^/"^ = V^. 
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5.2. Two-wheeled planar mobile robot 

Consider the motion of a two-wheeled planar mobile robot (or "two- wheeled carriage") which is 
able to move in the direction in which it points and, in addition, can spin about a vertical axis 



1 18, 20, Let P be the intersection point of the horizontal symmetry axis of the robot 

and the horizontal line connecting the centers of the two wheels. The position and orientation of 
the robot is determined, with respect to a fixed cartesian reference frame (with horizontal x- and 
y-axes) by {x,y,6) G SE{2), where E §^ is the heading angle and the coordinates (x,y) G 
locate the point P (see ^]). Let ipi, ^2 £ denote the rotation angles of the wheels which are 
assumed to be controlled independently and roll without slipping on the floor. The configuration 
space of this system is Q = x x SE{2). 

The Lagrangian function is the kinetic energy corresponding to the metric 

g = mdx ® dx + mdy ® dy + uiqI cos 9{dy ^ d9 + dO ^ dy) 

—mol sm9{dx ^ d6 + d6 ® dx) + JdO ® dO + J2d7pi dipi + J2d7p2 ® dijj2 , 

where m = niQ + 2mi, tuq is the mass of the robot without the wheels, J its moment of inertia 
with respect to the vertical axis, mi the mass of each wheel, Jw the axial moments of inertia of 
the wheels, and / the distance between the center of mass C of the robot and the point P. 

The constraints, induced by the conditions that there is no lateral sliding of the robot and that 
the motion of the wheels also consists of a rolling without sliding are: 

X sin — y cos 9 = , 
X cos 9 + y sin 9 + c9 + Ri^i = , 
X cos 9 + ij sin 9 — c9 + R'ip2 = , 

where R is the radius of the wheels and 2c the lateral length of the robot. The constraint distri- 
bution D is then spanned by 

d R f I d\ d R f Id 
cos 9— — hsmt'— — I — -— , — cosy-- — h sm 



5-01 2 \ dx dy c 89 J 5-02 2 \ dx dy c 89 

If we consider the action of G = SE{2) on Q 

* : GxQ — > Q 

((o, b, a), {ipi,tp2,x, y, 9)) 1 — > {ipi, tp2,a + x cos a — y sin a,b + x sin a -|- y cos a,a + 9) , 

we see that the constraint distribution V can be interpreted as the horizontal subspace of the 
principal connection 

(R R R R \ 

dx + — cos 9dipi + — cos 9dip2 + y(d9 + —dipi — —dip2) ] ei 
2 2 2c 2c J 

(R R R R \ 

dy + TT sm9d'tpi + — s'm9dip2 - x(d9 + —dipi - —d^2) ) 62 
2 2 2c 2c J 

R R 
+{d9 + — d0i - — dV'2)e3 , 
2c 2c 
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where {ei, 62, 63} is the canonical basis of the Lie algebra of G, with associated fundamental vector 
fields ^ ^ d d d 

(ei)Q = ^, (e2)Q = ^, (^3)Q = ^-y^ + x^. 

The curvature of 7 is 

^ = — — (sin0 ei — cos 9 e2)d'ipi A dip2 ■ 
The induced metric on M = Q/G is given by 

g = {J2 + m— + J—^)d'ilJi (g) di/ji 
4 4c^ 

^2 ^2 ^2 ^2 

+(m—; J—^)(dipi (g) ^-02 + d^/^2 ® dV'i) + (<^2 + m— + J— 5-)(iV'2 ® ^-02 • 

4 4c^ 4 4c^ 

The Lagrangian L* is the kinetic energy function induced by g. The gyroscopic 1-form otx here 
becomes 

molR^ ■ ■ ■ 

ax = ^^2 - '4^i)Widip2 - V'2'^V'i) ■ 

Then, the symplectic reduction (^) yields 

ixUJL* = dEi* +ax . 
On the other hand, the metric connection tensor K is given by 

K = — — idil^i dipi A d'ijj2 — dip2 ® dipi A dV'2) • 
4c^ 

It is easily seen that, in this case, the tensor field B is not skewsymmetric and so, by Proposition 
4.4, V^^"^ 7^ V^. In addition, the Christoffel symbols of the metric connection V are 

-pi'l -pi'2 TV- 

-pi'l _jy^ -pi>2 _ 

■p'01 _ _ ■p'/'2 _ _ 7^ 

tV"! _ r'/'2 _ TV- 

^ V'2t/'2 ~ -^2 , , J- V2l/'2 ~ ' 

where 

J2 — mc^ 



Ki = moR'^l- 



4c2(2Jc2 + 32R^){2J + mc^) 

K2 = R-'l- 



3 4Jc2 + (J2 +mc2)i?2^p 



'4c2(2Jc2 +j2^2)(2J + mc2 
Clearly, the torsion T does not vanish, and so V / V^. 
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6. Relation between both approaches 



The above examples show us the fohowing intriguing fact. In the case of the mobile robot with 
fixed orientation, the 1-form ax identically vanishes, and thus the reduced problem has no ex- 
ternal gyroscopic force in its unconstrained symplectic formulation. Consequently, since there is 
no potential, the solutions of the reduced system are geodesies of the Levi-Civita connection V^. 
Indeed, we verified that V = V^/^ = V^. However, in the case of the two-wheeled mobile robot we 
obtained Wx 7^ and V 7^ / V^/^ . Apparently, there exists a relation between the properties 
of the contorsions of the connections considered in Section 4 and the vanishing (or not) of the 
gyroscopic 1-form. 

Using the definition of ax in (0) , one can check that the following relation holds (see |jll| ) : for any 
Y € X{V), 

ax{Y) = yrr{Y){eL{X)) + OURiX, Y)) + eL{hr[X,^^r{Y)]) , 
where R is the tensor field of type (1,2) on V given by 

R=^[hr,hr] , 

with [ , ] denoting the Nijenhuis bracket of type (1,1) tensor fields. The relation between R and Jl^, 
the curvature tensor of the principal connection T, is given by R{U,V){v) = {Q'^ {Uv,Vu))tq{v) 
for any U,V € X(D) and v gV. 

In particular, if we take a horizontal vector field Y (^U, we deduce from the above that ax(Y) = 
9l{R{X,Y)). Herewith, the action of the gyroscopic 1-form 7xx on a vector field Z € X{T{Q/G)), 
evaluated at a point Wq G T{Q/G){= V), becomes 

a5^{Z){wg) = ax{Z''){vg) = {eL),^iRiX,Z^)ivq)) 

= (^l)., ((f^^(X.,,<))TQ(r',))=(eL)., ((l^^(TQ,X,„rQ^Z,^J)TQ(^;,)) , 

for an arbitrary Vq such that p{vq) = Wq, and where the last equality has been derived using 
In these expressions, Z^ is the horizontal lift of Z with respect to T. Recalling that the Poincare- 
Cartan 1-form 6^ and the Legendre mapping (induced by the given Lagrangian) FL : TQ — > T*Q 
are related by (^l)?,, (u) =< FL{vq),TQ^{u) >, for any u € T^^TQ, and taking into account that 
X is a SODE, we further obtain 

m{Z){wq) = <FL[vq),{W{TQ^X,^,rQ^Z':^))Q{q)> 



= gq(vq,{W{Vq,TQ^Z^^))Q{q)^ 
= g,{vq,{n^{Vq,{TQ/G^Z^,tq))Q{q)) . (13) 

Note that in the last expression, the horizontal lift of tq/q^Z^]^ is the one with respect to 7. 

An important observation is that (13) immediately shows that the gyroscopic 1-form oY is semi- 
basic with respect to the canonical fibration tq/q : T{Q/G) — > Q/G. Indeed, assume (tq/g')*oZ = 
0, then it readily follows that ax(Z) = 0. 
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Elaborating (|T^ ) a bit further, using the metric connection tensor K and the contorsion B intro- 
duced in Section 4.1, we find 

-a:^{Z){wq) = gg (vg,n{vg,{TQ/G^Z^-)'^))Q{q)^ 

This proves the next result, which was already implicit in the work of Koiller |20| ]. 

Proposition 6.1. An explicit relation between the gyroscopic 1-form and the contorsion tensor 
field B, defined in Section 4, is given by 

{ax)wg{u) = gq{B{Wq,Wq),TQlGu) , 

for all u G T^-T{Q/G)), Wg G T{Q/G). 

From this we can immediately deduce, taking into account Proposition 4.4: 
Corollary 6.2. The following statements are equivalent: 

(i) Wx vanishes identically; 
(a) B is skew-symmetric; 
(Hi) ^ y3_ 

Remark 6.3. One can think of simple examples in which B is skew-symmetric but nonzero. 
Consequently, if oY vanishes, this does not imply V = V^, although in such a case both connections 
do have the same geodesies (see Appendix). 

By means of Proposition 6.1, one can also recover the gyroscopic character of ax, established 
already in Proposition 3.4. For that purpose, let us define the following 2-form on T{Q/G): 

E{Y, Z){wg) = gq{B{wg,TQ/GY^^),rQ/G^Z^^) . (14) 

One readily verifies that H is indeed bilinear and by Proposition A. 3, H(y, Y) = 0. It is then easy 
to check that 

ax = ijc'^- 

In local coordinates q"" {a = 1, . . . , n — k) on M = Q/G, we have that 

H = ^ fB'^Jhcdq'' Adq\ «Y = ^ q'^fBl^c^M' ■ (15) 
a<6 a,e,c 

A careful calculation, very similar to the one performed for proving Proposition 6.1, reveals that 
for generalized Chaplygin systems of mechanical type, the 2-form S of Proposition 3.4 and the 
above 2-form H coincide. 
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7. Integrability aspects and the existence of an invariant measure 
7.1. K oilier 's question 

In [^], the author wonders whether there might exist an invariant measure for the reduced equa- 
tions of a generahzed Chaplygin system. In the following, we shah deal with this problem in some 
detail. 

The existence of a measure which is invariant under the flow of a given dynamical system is a 
strong property. Indeed, using an integrating factor it is possible to derive from it (locally) an 
integral of the motion. This fact plays an important role in discussions concerning the integrability 
of the system under consideration, as illustrated by the following theorem. 

Theorem 7.1.([^) Suppose that the system x = X(x), x £ N, with N an n-dimensional smooth 
manifold, admits an invariant measure and has n — 2 first integrals Fi, Fn~2- Suppose also that 
Fi, ...,Fn-2 are independent on the invariant set Nc = {x £ N : Fs{x) = Cg, 1 < s < n — 2}. Then: 

• the solutions of the differential equation lying on Nc can be found by quadratures. 

Moreover, if Lc is a compact connected component of the level set Nc and if X does not vanish on 
Lc, then 

• Lc is a smooth manifold diffeomorphic to a two-torus; 

• one can find angular coordinates (pi, (p2mod{2'K) on Lc in terms of which the differential 
equations take the simple form 

$((^i,(/?2)' ^"^ $((^l,(/?2)' 

where uji, lo2 are constant and ^ is a smooth positive function which is 27r -periodic in ipi, 

¥'2- 

By the Riesz representation theorem, we know that each volume form on an orientable manifold 
induces a unique measure on the Borel a-algebra [p. Therefore, with a view on tackling the 
integrability problem of generalized Chaplygin systems, it is worth looking for invariant volume 
forms under the flow of the reduced dynamics X. This is what we intend to do in the sequel. 

In this section, we again consider Chaplygin systems of mechanical type with a Lagrangian L = 
^g — V, where the metric g and the potential energy function V are G-invariant. The reduced 
equations of motion are (cf. (P)) 

ix^L* = dEi* + oY , (16) 

where L* = ^g — V, with g and V the metric and the potential function on M induced by, 
respectively, g and V. Remember that the energy El* is a constant of the motion. The local 
expression for the reduced dynamics takes the form (cf. (^)) 
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where F? are the Christoffel symbols of the Levi-Civita connection V^. 



The gyroscopic systems usually encountered in the mechanics literature |42, 43 1 differ in a crucial 
way from the ones we obtain through the reduction of a generalized Chaplygin system. In fact, 
the common situation in mechanics is that of a system, with configuration space P, described by 
an equation of the form 

where L : TP — M is a (regular) Lagrangian and where the gyroscopic force is represented by a 1- 
form a = i-p^TpH), with 11 a closed 2-form on P. These systems are then Hamiltonian with respect 
to the symplectic 2-form to = ujj^ — TpH, and thus they admit an invariant measure, determined 
by the volume form cj" = ujj^. 

In some sense, our reduced system ( |T6|) exhibits the opposite behaviour. Indeed, the 2-form H, 
defined by is semi-basic but in general it is not basic, i.e. it is not the pull-back of a 2-form 



on the base M = Q/G. This can be readily deduced from its local expression (cf. 15). Moreover, 
using (|l5|), the following property is easily proved. 



Proposition 7.2. The 2-form H is closed if and only if it is identically zero. 

Note, in passing, that a similar property also applies to the gyroscopic 1-form Sx- The semi-basic 
character of H ensures, however, that the 2-form lol* — H is still nondegenerate and, consequently, 
we have that the equation (|T^ ) can be rewritten in the form 

ixUJ = dEi* , (17) 

with Lo = ul* — H an almost symplectic form (i.e. a nondegenerate, but not necessarily closed 
2-form) . 



In IS^j, S.V. Stanchenko has studied Chaplygin systems of mechanical type with an abelian Lie 
group in terms of differential forms, in a way which shows many links to the approach described 
in Section 3. In our setting, his results can be generalized to the non-abelian case for any kind of 
generalized Chaplygin system. 



Let us assume, following |3^, that there exists a function F e C°°{T{Q/G)) such that 

dFAeL'=E. (18) 

Putting = expF, we have that 

diNuj) = d{NujL* - NE) = d{NLUL* -dN AOl-) = 0. 
Since (^) can still be written as 

"i-x/Ni^^) = dEi* , 
we deduce that Cxij^{Nuj) = 0. Consequently, 

= /:^/^(7V^)" = /:^iV"-ia;", 
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and we see that A^" ^cj" is an invariant volume form. This proves the following result. 

Theorem 7.3. (p^]) Condition ^TB^ is sufficient for the existence of an invariant measure for the 
reduced Chaplygin equations (|7^. 

Remark 7.4. Stanchenko observes that if F satisfies ([T8|), the semi-basic character of both 6^* 
and H imply that F is necessarily the pullback of a function on Q/G{= M). 

It turns out that condition (|18|) can be relaxed to some extent: it suffices to require the almost 
symplectic 2-form uj to be globally conformal symplectic, that is, that there exists a function 
F G C7°°(T(Q/G)) such that 

dF Auj = -duj . (19) 

Theorem 7.3 still holds in this case, with ( |l9|) replacing The previous remark also remains 

valid: the function F is necessarily the pullback of a function on Q/G. Note that ( p!^ ) obviously 
implies (|l9|). 

However, even the weaker condition (|T9| ) is not necessary in general for the existence of an invariant 
volume form on T{Q/G). To derive a necessary condition, let us suppose that fj, is an invariant 
volume form for the dynamics X on T(Q/G). We then necessarily have that fj, = kco'^, for some 
nowhere vanishing function k E G°°{T{Q /G)). Restricting ourselves to a connected component of 
Q/G need be, we may always assume k is strictly positive. It follows that 

= X{k)uj'^ + nkCxi^ = X{k)uj'^ - nkixdE Alo'^-^ . 

The 2n-form i^dE A uj^~^ determines a function h e G°°{T{Q/G)) by 

ijfdHAu;"-^ = ^w". (20) 

Therefore, we have 

X{k) = kh or, equivalently, X(lnA;) = h. (21) 

This essentially yields the same characterization as the one derived in [^]. Now, conversely, 
assume there exists a function k satisfying (21), with h defined by (^). Going through the above 
computations in reverse order then shows that the 2n-form /cu;" is an invariant volume form of X. 
We may therefore conclude that the existence of a globally defined function k for which (^Tj) holds 
is not only a necessary but also a sufficient condition for the existence of an invariant volume form. 
It is interesting to note that in |3^, Stanchenko has proved that in case the reduced Lagrangian 
is of kinetic energy type, L* = ^g, and if there exists a solution k of (^) which is basic, i.e. which 
is the pullback of a function on the base space Q/G, then the volume form = kiv'^ remains an 
invariant of the reduced dynamics if a potential energy function V € G°°{Q/G) is included in the 
Lagrangian L* (coming from a G-invariant potential added to the given Chaplygin system). 

Obviously, however, ( ]2l| ) is not a very handy criterium to deal with in practice. We will now see 
that, at least for a subclass of generalized Chaplygin systems, it may be replaced by a more easily 
manegeable condition. 
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From (po|), we can deduce a local expression for h. After some computations, we get 



a,b 



and, using (15), this further becomes 



f,b c,e 

Note that S*dh is the pullback of a basic 1-form, i.e. S*dh = TqjqI^^ where the local expression for 
fj reads 



P = he{q)dq- = Y.{Bl, + Bl,)dq' 



with he = dh/dcf. Let us assume now that there exists a basic function k for which ( piD holds. 
We then have that S*{dX{\nk)) = dihik). Therefore, taking the differential of both hand-sides of 
(|2l|) and applying S* to the resulting equality, we obtain 

d{\nk) = (3, 



or, in local coordinates. 



d{lnk) 



We thus see that, if (|2l| ) admits a solution k which is basic, then the 1-form /? is exact. 

It turns out that, for systems for which V = (i.e. L* is a pure kinetic energy Lagrangian), the 
previous result even has a more definitive character. Indeed, let V = and suppose there exists a 
function k E C°°{T{Q/G)) (not necessarily basic) satisfying (21). Then, we have that 



S*dX{\iik) = S*dh. 



In local coordinates, this becomes 



dX[\nk) 
dq^ 



he{q) , e = 1, ...,n. 



,,a(lnfc) ^ ^^,a(lnfc) 



But X(ln A:) = n^'r.illll^ + j^a ^y-yj ^^g^g ^ -(g"-^ab + fl^flT? ), and so we have 



dX{\ak) d{\xik) .^92 (In A;) dX"" d{\nk) -^d'^{\nk) 

+ q „ „ ^ ■„ H TTTZ TTTZ 1" X 



dq'^ dq'^ dq°-dq^ dq^ dq°- 

In points Og of the zero section of T{Q/G) this reduces to 

dX{lnk) d{lnk) 



QqaQqe 



dq^ 
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If we now define the basic function k = k o s, where s : Q/G — > T{Q/G), q i — > (g, 0) determines 
the zero section, we derive from the above that 

=^W = ^^(,.0)=M.).e=l....,n. 

and, hence, it follows again that the 1-form /? is exact. 

Conversely, if the 1-form /3 is exact, say (3 = df for some function / € C°°{Q/G), and putting 
k = exp(TQ/G*/), one easily verifies that k satisfies (pl|). This obviously also holds in the presence 
of a potential (i.e. V ^ 0). 

Summarizing the above, we have proved the following interesting result. 

Theorem 7.5. For a generalized Chaplygin system with a Lagrangian of kinetic energy type, i.e. 
L = ^g, there exists an invariant volume form for the reduced dynamics X on T(Q/G) iff the basic 
1-form P, defined by S*dh = Tq^^(3 (with h given by ^dj)), is exact. The 'if part also holds if L 
is of the form L = ^g — V , with V a G-invariant potential. 

Therefore, if we manage to find a particular example of a system with a kinetic energy type 
Lagrangian for which (3 is not exact, we shall have proved that the answer to Koiller's question 
about the existence of an invariant volume form for all generalized Chaplygin systems, is negative. 
In particular, for such a counter example it suffices to show that the corresponding [3 is not closed. 

7.2. A counter example 

Let us consider the following modified version of the classical example of the nonholonomic free 
particle ^j. 

Consider a particle moving in space, so Q = M'^, subject to the non-holonomic constraint 

(j) = z — yxx . 

The Lagrangian function is the kinetic energy corresponding to the standard metric g = dx®dx + 
dy ® dy + dz dz. Therefore, 

L = \{x' + y^ + i^) . 

The constraint submanifold is defined by the distribution 

^ d d d 

ox oz oy 

Consider the Lie group G = M with its trivial action by translations on Q: 

$ : G X Q — > Q 

{s,{x,y,z)) I — > {x,y,z + s). 

Note that V is the horizontal subspace of a connection 7 on the principal fiber bundle Q — > Q/G, 
where 7 = dz — yxdx. Therefore, this system belongs to the class of generalized Chaplygin systems. 
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The curvature of 7 is given by 

^l'^ = xdx A dy . 

The induced metric g on Q/G = is 

^ = (1 + x^y^)dx ® dx + dy dy . 

The metric connection tensor K is determined by 

^(959 2 Ki ^ ^ ^ ^ 
dx^ dx^ dy ^ ^ ' dy' dx' dy 

Then, the contorsion of the affine connection here reads 

2 , T "9 x^y , , d 

B = X ydx ® dx (g) — Tr-^dx ^ dy (g) — . 

dy 1 + x^y^ dx 

Finally, the 1-form /? associated to the reduced Chaplygin system is given by 

^ = — T^^y^ 

1 + x^y^ 

which is clearly not closed. Hence, according to Theorem 7.5, there is no invariant volume form 
for the reduced dynamics. 

Note that in this example the distribution D has 'length' 1 at all points of not belonging to the 
plane x = 0, since T> + [V, V] spans the full tangent space at all points for which x 7^ 0. We leave 
it as a challenge for the reader to find a more generic example (on M'^) for which the constraint 
distribution has length 1 everywhere. 

Remark 7.6. The classical model for the nonholonomic free particle corresponds to the constraint 
(f) = z — yx. Here the constraint distribution does have length 1 everywhere, however, after 

y 

performing the appropriate computations, we obtain (3 = ndy, which is clearly exact: (3 = 

1 + r 

df, with f = —In I rr I . Then, k = — , and the invariant measure defined by kuj" leads, 

2 \l + yy ^/TT7 
using Euler's integrating factor technique, to the constant of the motion 



ip = Xy/l + y2 

which was also obtained by different methods in ^, ^, p^ . 

Appendix: Metric connections 

In this section, we want to collect some simple facts about metric connections that have been useful 
for the formulation of generalized Chaplygin systems. 

Let Q be a n-dimensional differentiable manifold, the configuration space with Riemannian metric 
g. We denote by the Levi-Civita or Riemannian connection associated to the metric g. 
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Definition A.l. An affine connection V is called metric with respect to g ifVg = 0, that is, 

Z{g{X, Y)) = g{VzX, Y) + g{X, VzY) , 

for all X,Y,Z e X{Q). 

Let V be a metric connection with respect to g. In the fohowing proposition, we prove that V is 
determined by its torsion T. 

Proposition A. 2. Let T be a skewsymmetric (1,2) tensor on Q. Then there exists a unique 
metric connection V whose torsion is precisely T. 

Proof. Let us suppose that there exists such metric connection V. Then we have that 

Z{g{X,Y)) = g{VzX,Y) + g{X,VzY), 
X{g{Z,Y)) = g{VxZ,Y)+g{Z,VxY), 
Y{g{X,Z)) = g{VYX,Z)+g{X,VYZ), 

for all X,Y,Z e X{Q). Now 

Zig{X, Y)) + X{g{Z, Y)) - Y{g{X, Z)) 
= g{VxZ + VzX, Y) + g{VzY - VyZ, X) + g{V xY - VyX, Z) 
= g{2VxZ + T{Z, X) + [Z, X],Y) + g{T{Z, Y) + [Z, F], X) + g{T{X, Y) + [X, Y], Z) 
= 2g{VxZ, Y) + g{T{Z, X) + [Z, X],Y) + g{T{Z, Y) + [Z, Y],X)+ g{T{X, Y) + [X, Y],Z). 

Consequently, the connection V is uniquely determined by the formula 

g{VxZ, Y) = 5(V^Z, Y) - ^ ig{Y, T(X, Z)) + g{X, T{Z, Y)) + giZ, T(X, Y))) . (1) 

\qed\ 

This proposition implies that the Christoffel symbols T^c metric connection V in a local 

chart (q^) are given by 

^BC = Tfic - 2^^^ {9kmTbc + 9bmTck + 9cmTbk) > 
where are the Christoffel symbols of the connection and T = T^^dq^ ^ dq^ ® 

Another way to characterize metric connections is the following. Given any affine connection on 
Q, we know that 

VxY = V'xY + SiX,Y), 
where 5 is a (1,2) tensor field, called the contorsion. If V is a metric connection, then 

Z{g{X,Y)) = g(yzX,Y)+g{X,VzY) 

= 5(V|X + S{Z, X), Y) + g{X, V|y + S{Z, Y)) 

= Z{g{X, Y)) + g{SiZ, X),Y) + g{X, S{Z, Y)) , 
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which imphes that g{S{Z, X), Y) + g{X, S{Z, Y)) = 0. Then we have proved the following 
Proposition A. 3. V is a metric connection if and only if 

g{S{Z,X),X) =0 yX,Z e X{Q) . (2) 

As a consequence of the two characterizations we have obtained for metric connections, we can 
establish the next result. 

Corollary A. 4. There is a one-to-one correspondence between (1,2) tensors S verifying and 
skewsymmetric (1,2) tensors T. This correspondence is given by 

S — >T, 

where T{X, Y) = S{X, Y) - S{Y, X) and 

T — > S, 

where g{S{X, Z), Y) = {{g{Y, T{X, Z)) + g{X, T{Y, Z)) + g{Z, T{X, Y))). 
The equations for the geodesies of a metric connection are 

V,(,)c(t) = ^ V^(,)C(t) = -S{c{t), d{t)) , 

or, in local coordinates 

■■A I -pA ^B-C _ \ " „AK / rpM I „ rpM \ -B-C 

Q + ^ BcQ 1 -2^9 [gBivUcK + 9cmJ^bk) Q Q ' 

B<C 

for each A = 1, n. 

Finally, it is very important to note that the metric connections preserve obviously the kinetic 
energy K of the metric g, that is, if c{t) is a geodesic of V, we have that 

^ Q5(c(t),c(t)^ =5(Ve(t)c(t),c(t)) = 0. 
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